ON THE EXISTENCE OF SOLITARY TRAVELING WAVES FOR 
GENERALIZED HERTZIAN CHAINS 



ATANAS STEFANOV AND PANAYOTIS KEVREKIDIS 

Abstract. We consider the question of existence of "bell-shaped" (i.e. non-increasing 
I for a; > and non-decreasing for a; < 0) traveling waves for the strain variable of the 

' generalized Hertzian model describing, in the special case of a p = 3/2 exponent, the 

dynamics of a granular chain. The proof of existence of such waves is based on the 
English and Pego [Proceedings of the AMS 133, 1763 (2005)] formulation of the problem. 
More specifically, we construct an appropriate energy functional, for which we show that 
the constrained minimization problem over bell-shaped entries has a solution. We also 
provide an alternative proof of the Friesecke-Wattis result [Comm. Math. Phys 161, 394 
(1994)], by using the same approach (but where the minimization is not constrained over 
bell-shaped curves). We briefly discuss and illustrate numerically the implications on the 
doubly exponential decay properties of the waves, as well as touch upon the modifications 
of these properties in the presence of a finite precompression force in the model. 
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1. Introduction 
V ' 

, Localized modes on nonlinear lattices have been a topic of wide theoretical and exper- 

^ I imental investigation in a wide range of areas over the past two decades. This can be 
^ ■ seen, e.g., in the recent general review [1], as well as inferred from the topical reviews in 
■ nonlinear optics [2j, atomic physics [3] and biophysics [4] where relevant discussions have 
Q ! been given of the theory and corresponding applications. 

One of the areas in which the theoretical analysis has been especially successful in 
describing experimental data and providing insights has been that of granular crystals [S] . 
These consist of closely-packed chains of elastically interacting particles, typically according 
^ ! to the so-called Hertz contact law. The broad interest in this area has emerged due to the 
wealth of available material types/sizes (for which the Hertzian interactions are applicable) 
and the ability to tune the dynamic response of the crystals to encompass linear, weakly 
nonlinear, and strongly nonlinear regimes O |6l [71 E] . This type of flexibility renders these 
crystals perfect candidates for many engeenering applications, including shock and energy 
absorbing layers [HI [101 El [12], actuating devices [13\, and sound scramblers [HI [TJ]. It 
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should also be noted that another aspect of such systems that is of particular appeal is 
their potential (and controllable) heterogeneity which gives rise to the potential not only 
for modified solitary wave excitations [16], but also for discrete breather ones fTT]. 

Another motivation for looking at waves in such lattices stems from FPU type problems 
[m [19]. In the prototypical FPU context, it has been rigorously proved that traveling 
waves exist which can be controllably approximated (in the appropriate weakly supersonic 
limit) by solitary waves of the Korteweg-de Vries equation [20]. However, in more strongly 
nonlinear regimes, compact-like excitations have been argued to exist [5l [6l [8] (see also pT] 
for breather type excitations) and have even been computed numerically through iterative 
schemes [221 E3], but have not been rigorously proved to exist in the general case. In the 
work of [21], the special Hertzian case was adapted appropriately to fit the assumptions 
of the variational-methods' based proof of the traveling wave existence theorem of f25| in 
order to establish these solutions. However the proof does not give information on the 
wave profile. 

Our aim herein is to provide a reformulation and illustration of existence of "bell-shaped" 
traveling waves in generalized Hertzian lattices. Our work is based on the iterative schemes 
that have been previously presented in [23l [22] for the computation of traveling waves in 
such chains of the form: 

(1) = [Vn-l - VnY^ - [Vn - V„+l]^. 

Here f„ denotes the displacement of the n-th bead from its equilibrium position. The 
special case of Hertzian contacts is for p = 3/2, but we consider here the general case of 
nonlinear interactions with p > 1. Notice that the "+" subscript in the equations indicates 
that that the quantity in the bracket is only evaluated if positive, while it is set to 0, 
if negative (refiecting in the latter case the absence of contact between the beads). The 
construction of the traveling waves and the derivation of their monotonicity properties will 
be based on the strain variant of the equation for m„ = t'„„i — f„ such that: 

(2) Un = [Un+l]+ - 2[Un]+ + [Un-l] + , 

Our presentation will proceed as follows. In section 2, we will give a preliminary math- 
ematical formulation to the problem, briefiy illustrate its numerical solution and some of 
its consequences. Then, we will proceed in section 3 to state and prove our main result. 
Some technical aspects of the problem will be relegated to the appendices of section 4. 

2. Preliminaries and Numerical Results 

When seeking traveling wave solutions of the form m„ = u{x) = u{n — ct), we are led to 
the advance-delay equation (setting c = 1) 

(3) u"{x) = uP{x + 1) - 2uP{x) + uP{x - 1). xe 

where m is a smooth and positive function, with the desired monotonicity involving decay 
in (0, oo) and increase in (— oo,0). 



SOLITARY WAVES FOR HERTZIAN CHAINS 3 

2.1. Fourier transform and Sobolev spaces. We introduce tlie Fourier transform and 
its inverse via 



/oo 
-oo 

/oo 
-oo 



As is well-known, the second derivative operator has a simple representation via the 
Fourier transform, namely 



For every s > 0, we may define (— = {^'^'^^'^Y fiO the Sobolev spaces W^'P 
via 

ii/ik- = ii/iiL^ + ii(-Ar/viu.. 

We will also consider the operator 

Ad^scf{x) = fix + 1) - 2/(x) + fix - 1) 

on the space of L'^{'R^) functions. Using Fourier transform, we may write 

/oo pOO 
/(0(e^"'^ + e"^"'^ - 2)e2™«c/e = sin2(7rO/(e)e2™«cie 
oo J — oo 

In other words, ^disc is given by the symbol — 4sin^(7r^), that is 

(4) A;;~J(0 = -4sin2(7rO/(0- 

2.2. The English-Pego formulation. We may rewrite the equation (j3]) in the form 

(5) u"ix) = /\dUu'']ix) 

Taking Fourier transform on both sides of this (and using (jlj)), allows us to write 
-47r2e2M(0 = -4sin2(7rO^>(0 or 

Equivalently, taking A : A(^) = ^^^jp''' , 

/oo 
Kix - y)u^iy)dy =: M[u% 
-oo 

In other words, we have introduced the convolution operator M. with kernel A(-). It is 
easy to compute that A(x) = (1 — or 



(6) m) = '-^um 



A(x) 



1 — |x| |x| < 1, 
\x\ > 1. 



Note that we have the following formula for the convolution A * / 
(8) Mf = A* fix) = -\x- y\)fiy)dy. 

Jx-l 
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2.3. Numerical computations and other consequences of the English-Pego for- 
mulation. For reasons of completeness and in order to appreciate tlie form of (suitably 
normalized) solutions of Eq. ([7]), in Fig. [H we used this equation as a numerical scheme 
and proceed to iterate it until convergence. The figure illustrates the converged profile 
of the solution and its corresponding momentum 0j = —apx (for c = 1). The results of 
these computations are shown for different values of p (in order to yield a sense of the 
p— dependence of the solution, namely for p = 3/2 (the Hertzian case), p = 2 and p = 3 
(the FPU-motivated cases, in that they are the purely nonlinear analogs of a- and /3-FPU 
respectively) and finally p = 10 (as a large-p case representative). The figure shows the 
solutions' profile and corresponding momenta, as well as the semi-logarithmic form of the 
profile, so as to clearly illustrate the doubly exponential nature of the decay (see below). 
Notice that as p increases, the decay becomes increasingly steeper. 

To corroborate the exact nature of such traveling wave solutions, once the solution 
was obtained, then the "lattice ordinates" of both the solution and its time derivative 
were extracted and inserted as initial conditions for the dynamical evolution of Eq. (|2]). 
The results of the relevant time integration (using an explicit fourth-order Runge-Kutta 
scheme) are shown in Fig. [2l It can be straightforwardly observed that excellent agreement 
is obtained with the expectation of a genuinely traveling (without radiation) solution with 
a speed of c = 1, so that its center of mass moves according to x = ct (the solid line in 
the figure). This confirms the usefulness of the method (independently of the nonlinearity 
exponent p, as long as p > 1) in producing accurate traveling solutions for this dynamical 
system. 

If the convergence to such a nontrivial profile is established (as we will establish it 
in section 3 with the proper monotonicity properties based on our modified variational 
formulation), there is an important immediate conclusion about the decay properties of 
such a profile. In particular, 

(9) u{x + l) = j K{y)u'\x + l-y)dy <vF{x) ^ u{x + n) <u{xY" , 

Hence, as was originally discussed in [26] and then more rigorously considered in [22] (see 
also P3|), the solutions of Fig. [T] feature a doubly exponential decay. This very fast decay 
(and nearly compact shape) of the pulses can be clearly discerned in the semi-logarithmic 
plots of the figure. 

As a slight aside to the present considerations, we should mention that a physically 
relevant variant of the problem consists of the presence of a finite precompression force Fq 
at the end of the chain [3 [6]. In that case, the model of interest becomes (in the strain 
formulation and with Fq = 5q) 

(10) iln = [5o + M„+i]+ - 2[5o + + [5o + 

The case of 5o = constitutes the so-called sonic vacuum [5J, while that of finite 6q 
features a finite speed of sound (and allows the existence and propagation of linear spectrum 
excitations). It is worthy then to notice that for 5o 7^ 0, the above decay estimate is 
modified as: 

(11) u{x + 1) ~ 6q~^u{x) =^ u{x) ~ (5o~"^ exp(n logr(xo)) 
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Figure 1 . Each one of the panels illustrates the numerically exact (up to a 
prescribed tolerance which we set here to lO"*^) solution profile of the iterative 
scheme, as renormalized for use in Eq. (|2]). The solid (blue) line illustrates 
the spatial form of the solution and the dashed (red) line the corresponding 
momentum (for speed c = 1). The circles and stars denote respectively 
the ordinates of the lattice nodes (extracted for use in Eq. ([2]). The inset 
illustrates the profile in a semilog to highlight the doubly exponential nature 
of the decay (notice also the steepening as p increases). The top left panel 
is for p = 3/2, the top right for p = 2, the bottom left for p = 3 and finally 
the bottom right for p = 10. 

Namely, the solutions are no longer doubly exponentially localized but rather feature an ex- 
ponential tail (and are progressively closer to regular solitary waves). This can be thought 
of as a "compacton to soliton" transition that is worth exploring further (although the 
case of (5o 7^ will not be considered further herein). 

2.4. Some facts and definitions regarding distributions. We now turn to several 
definitions, which will be useful in the sequel. Our space of test functions will be the 
following. For C - an open set, let VCV) = C^{V) be the set of all C°° functions 
with compact support, contained inside V. We equip this with the usual topology of a 
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Figure 2. Result of direct integration of Eq. ([2]), with r„(0) and r„(0) as 
seeded from the iteration scheme's convergent profile. The solid line in each 
case illustrates the trajectory of x = ct (for c = 1) to which the solutions 
correspond. One can notice for all values of p (p = 3/2: top left; p = 2: top 
right; p = 3 bottom left and p = 10 bottom right) the agreement with the 
expectation of a genuinely traveling (non-radiating) waveform of c = 1. 



Frechet space, generated by a family of seminorms PnU) = sup^gy^ Ylia=i \^xfi^)\i where 
is some fixed nested family of compact sets, so that U„\4 = V . The distributions 
over this space of functions, which we denote by ViV), is its dual space, namely all 
continuous linear functionals over ViV). The derivatives of such distributions are defined 
in the usual way {h'^ip) := —{h,tp'). One may also define a convolution of a distribution 
h G V{lV^) with a given C°° function x, by {h * x, ip) := {h, ip * x{~'))- 

We say that two distributions are in the relation hi < h2 in V{V) sense, if for all 
ip e V{V),ip > 0, we have {hi,ip) < {h2,ip). In particular. 

Definition 1. We say that a distribution h is non-increasing (non-decreasing) over a set 
V, ifh'<0 (h' >0) m V'{V) sense. 
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Of course, if h happens to have a locally integrable derivative h' on an interval, then the 
notion of non- increasing function coincides with the standard (pointwise) notion by the 
fundamental theorem of calculus. More generally, we have the following 

Lemma 1. Suppose that h is a locally integrable function in {a,b) and it satisfies h' > 
(h' < in V'{a,b) sense. Then, h is almost everywhere (a.e.) non- decreasing (non- 
increasing, respectively) function on (a, b). That is, for almost all pairs x < y, h{x) < h{y) 
(h{x) > h{y) respectively). 

Proof. It is well-known by the Lebesgue differentiation theorem that for a locally integrable 
function h, one has h{x) — lim5_^o | J^^^ h{y)dy a.e. All such points x are called Lebesgue 
points for h. Denote this full measure set by L. We will show that for all a < Ci < C2 < 6 : 
Ci, C2 G L, we have h{ci) < h{c2). Indeed, let (5 > be so that S < min(c2 — Ci, Ci — a, 6 — C2). 
Define a function tjj, 

X < Ci — 5 

X — {ci — S) ci — 5 < X < ci 

1p = < S Ci < X < 02 

C2 + S — X C2 < X < C2 + S 

C2 + 5 <x. 

Clearly ^ > is not smooth, but is continuous and it can be approximated well by test 
functions. Moreover ip' = 1 on (ci — S, Ci), ip'i^x) = —1 on (c2, C2 + S) and zero otherwise. 
Since < {h',ip) — —{h,ip') (and these are well-defined quantities), we obtain 

{h,^l^')= r h{y)dy- r^\{y)dy 

J Ci—S J C2 

This is true for all 6 > 0, which are sufficiently small. Thus, dividing by 5 > and taking 
limit as (5 — )■ 0+ (and taking into account that both ci, C2 are Lebesgue points), we conclude 
that h{ci) < h{c2). □ 

We find the following trick useful, which allows us to reduce non-increasing/ non- decreasing 
distrubutions to non-increasing/non-decreasing functions. 

More precisely, let us fix a positive even function $, so that supp $ C (—1, 1) and 
J ^{x)dx = 1. Let ^s{x) ■— |'^*(f ) and define hs — h* ^s- The following lemma has a 
standard proof. 

Lemma 2. Let h be a non-increasing (non-decreasing) distribution. Then for every S > 
0, hs is a function, which is non-increasing (non- decreasing respectively). Moreover 
hs ^ h in the sense of distributions. 

Next 

Definition 2. We say that a distribution h G V'(Ti^) is bell-shaped, if h is non- decreasing 
in P'(— 00, 0) and h is non-increasing in V'{0, 00). 

In the sequel, we need the following technical result. 

Lemma 3. Suppose that f is an even distribution, so that f is non-increasing in V'^O, 00) 
and non- decreasing in V{—oo,0). Then A * f is non-increasing in (0, 00) and non- 
decreasing in (—00,0). 
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Assume that g is non-increasing in T''(a, oo) for some a G R^. Then A * g{x) is non- 
increasing in {a + 1, oo). 

Assume that h is non- decreasing in (—00, a) . Then, A*h is non- decreasing in (—00, a — 

l). 

Proof. By Lemma [21 it suffices to consider functions instead of ditributions with the said 
properties. We have the following computation for the derivative of the function A* g 

(12) (A*gy{x)= / giy)dy- / giy)dy 

J X J x — 1 

which follows by differentiating (|8]). 

It is immediate from ( 1T2|) that the claims for the functions A * g and A * h hold true. 
Regarding A * /, it is clear that A * / is non- increasing in (1, 00) and non-decreasing in 
(—00, —1). Thus, we need to show that for x G (0, 1), (A * g)'{x) < and for x G (—1, 0), 
(A * g)'{x) > 0. We only verify this for x G (0, 1), since the other inequality follows in a 
similar manner. Indeed, using f{x) = /(— x), 

(A*/)'(x)= r'f{y)dy- r f{y)dy= f^' f {y)dy - { f f {y)dy + f ^ f{y)dy). 

Jx Jx-1 Jx JO Jo 

Since / is non- increasing in (0, 00), 

ex px /•2x 

f{y)dy > f{y + x)dy = / fiy)dy, 

Jo Jx 

' ^ f{y)dy> f ^ f{y + 2x)dy= f^' f{y)dy 

-'0 J2x 

Going back to the expression for (A * /)'(x), this implies (A * fYix) < 0, which was the 
claim. □ 

We will also need the following multiplier in our considerations 

It is actually easy to see that since = ~|^) we have also the representation 

(13) Qf{x)= r'^'fiy)dy. 

Jx-1/2 

Noted that by our definition of the operator M., we have Ai = Q^. 

3. Statement and Proof of the Main Result 

Theorem 1. The equation (jS]) has a positive solution u, so that 

• u is even, 

• u is bell-shaped, 

• M G H°°(R^) 
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Notice that the fundamental contribution of our work to the setting of generahzed 
Hertzian lattices is the characterization of the (monotonically decaying on each side) na- 
ture of the traveling waves. We first explain the possible approaches to this problem. One 
approach is to use the form (IHl), to show that the map — )• A * 0^ has a fixed point. We 
shall pursue a different route, more in line with the work of Friesecke and Wattis [25] . 
Namely, we shall consider a constrained minimization problem associated with (|3]). We 
should note here, that the existence result of Friesecke- Wattis is based on the equivalent 
via the change of variables u{x) = z{x — 1) — z{x) (as discussed in section 1) formulation 

(14) z"{x) = [z{x - 1) - z{x)]l - [z{x) - z{x + 

We will in fact consider a different representation of (j5]). To that end, introduce a 
positive function w : w^^^ = u, whence ([5]) reduce£| to 



in^ (vr, 



This is easily seen to be equivalent to w^/p{^) = ^^^375^'u)(0- Undoing the Fourier transform 
yields 

(15) w^/P = Mw = Q'^w = A*w. 

Thus, we need to find a solution w to f|T5|) . which is as stated in Theorem [H 

Let e > be chosen appropriately small momentarily. Let g = 1-|-^ G (1, 2) and consider 
the following constrained optimization problem 



(16) 



Je{v) = /f -1 \Qv{x)\'^dx — max 
subject to I{v) = f^iv'^{x)dx = 1, 
V > 0,v — bell-shaped 



Solving (fT6l) will eventually lead us to a solution of (fT5|) . 

3.1. Constructing a maximizer for (ITB]) . Let us first show that the expression Js{v) is 
bounded from above, if v satisfies the constraints. Indeed, we have by Sobolev embedding 

Je{v)<\\Qv\\l,<C,\\Qv\\%,,^.,,,,^. 

By Gagliardo-Nirenberg inequality, we further bound 

IIQ^IIv.v,-v.,. < \\Qv\\XZ'^'\\Qv\\%'-'/'. 

But by the definition of Q (see ( |T3l) ). we have 

IIQ^IUm = \mQv]\\L^ = H- + 1/2) - vi- - 1/2)|U, < 2\\v\\l, = 2, 

r.Z+l/2 



(/ viy)dyydx<\\v\\% = l. 

OO Jx-1/2 



Thus, Je{v) < 2Cq and hence J^^^ = sup^, J(f ) < 2Cq, where Cg is the square of the 
Sobolev embedding constant W^^''-^/^'i(Tl^) ^ L'^(B}). 



^Note that this is a good transformation, since we are interested in positive solutions of 
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We will now select Eq so small that J^^^ = sup^ J{v) (which we have just shown exists) 
is positive for all < e < Eq. Indeed, take vo{x) = Cqe~^ , so that c'^fT^e~'^^ dx = 1. 
Thus, the function vo satisfies the constraints and therefore 

Jr"" > Avo) = cl I / e-^y"dy\^dx. 

Clearly, as e — )• 0, we have that the right hand side converges to 

I ^~'^''^^dy?'dx > 0. Thus, there exists Eq, so that for all e e {0,Eo), J^""" > 0. 

In fact, we can select the Eo, so that 

„2 poo px+l/2 

(17) J,'""" > ^ / I / e-^y^dy\''dx, 

2 J-oo Jx-l/2 

whenever e < Eq. For the rest of this section, fix e < Eq. Construct a maximizing sequence 
v"- e L^(R^) : ||w"||l9 = 1, so that Jeiv"^) — )■ J™"""". Namely, we take t;" satisfying the 
constraints, so that Jeiv'^) > J^^^ — ^. We only consider n large enough, so that ^ < J^^^, 
in which case Je(f") > 0. 

By the compactness of the unit ball of L'^ in the weak topology, we may take a weak 
L'^ limit V = lim„f„. Clearly, weak limits preserve the property that v is even and that v 
is a bell-shaped function. We now need to show that v satisfies the constraint \\v\\l<i = 1, 
which is non-trivial since norms are in general only lower semicontinuous with respect to 
weak limits (and hence, we can only guarantee ||f < !)• 

We show now that, there exists a subsequence {nk}k so that 

\Qv'"'ix)\^dx= / \Qvix)\^dx 

Indeed, we check that 

WQv^w^. = mQv^]\\L, + WQv'^Wl. = \\v\ + 1/2) - t;"(- - 1/2)|U, + WQv'^h. 

Thus, by the compactness of the embedding W^''^(R}) d L'^(—e~^,e~^) we conclude that 
there is a subsequence {nk}k and z G L^, so that HQf"*" — z\\l2 — )• 0. By uniqueness of 
weak limits, (note that Qv^'' Qv), it follows that z = Qv and hence f|T8l) . 
Clearly now 

\Qv'"=ix)\^dx) = ( / \Qv{x)\^dx) = Je{v). 

-£-1 

On the other hand, ||f ||l9 < 1, by the lower semi continuity of || ■ with respect to 
weak limits, but clearly f 7^ 0, since Je{v) = J^^^ > 0. We will show that in fact 
\\v\\iq = 1. Indeed, assume the opposite < p = \\v\\li < 1 and consider the function 
v/p : ||f/p||L9 = 1. Observe that 

Thus, ||f IIl"? = 1 (otherwise, we get a contradiction with the constrained maximization 
problem (fT6|) ). This implies that Je{v) = J™*^^, otherwise, we get a contradiction with the 
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definition of J™^^. Tlius, we liave sliown tliat tlie weak limit v is indeed a maximizer for 



3.2. Euler-Lagrange equations for the maximizer of f|T6|) . Before we proceed to tfie 
actual proof, it is relevant to mention a few words about our strategy. We will have no 
essential difficulties in deriving the Euler-Lagrange equation (see (1201) below) on the selB 
Q ^ {x : v'{x) < 0} using the standard calculus of variations arguments. The reason is that 
on the compact subsets of Q, the maximizer v is already a strictly decreasing function and 
therefore, for each fixed z, with support inside fl and regardless of its increasing/decreasing 
behavior, there will be A = A(2;) > 0, so that v + Xz will be an acceptable function for the 
maximization problem (fT6!l . 

On the other hand, we will have issues deriving appropriate equations on \ Q, the 
reason being that on every non-trivial interval of the set, we will have v = const, and hence 
the perturbations z must be increasing in this interval, in order for v + Xz, A > to be an 
acceptable function for the maximization problem (1161) . 

We will show that the set \ f2 consists of isolated points. To that end, under the 
assumption that there are non-trivial intervals inside R^ \ Q, we derive an Euler-Lagrange 
equation for v on such intervals (see (1291) below), which in turn will imply that the maxi- 
mizer V is trivial, a contradiction. 

First, we start with a technical issue. Since v > and f J, in (0, oo), there is a = inf{x > 
: v{x) = 0} < oo. If cr < oo, we henceforth restrict our attention to (— cr, cr). Clearly 
f^{y) = for all y > a. 

Consider perturbations of v of the form v + Xz, where A > and z G C^(R^), so that 
supp z C {—£~^ + 1, — 1) n (— cr, cr). We have 



We now define the set Vt. Roughly speaking, we would like to define 

Vl = {x : v'{x) < 0}. This is however impossible, since v' is merely a distribution. 

Instead, we define to be the maximal open subset of (0, oo), so that for every com- 
pact subinterval of it [c,d] C Q^, there is 5c,d > with v'{x) < —Sc^d in V^c^d) sense. 
Equivalently, we may define 



(USD. 



J — oo 




Je{v + A^;) 



{Qiy + Xz), Q{v + A2;))l2(_£-i^£ 
J,{v) + 2X{Q^v,z) + 0{X^). 



{xo > : 3ro, 6 > 0, {v', ^) < -5 / ^ : V^^ > G I?(xo - ro, Xq + ro)} 



Define Vt' = {-x : x G f^+j and finally = fi" U 1]+. 



■for the precise definition of see below 
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3.2.1. Euler- Lagrange on Vt. Due to our requirement for bell-shaped test functions in f[T6|) . 
we need to impose extra restrictions on the function z. To that end, note that VL is an 
open set and fix an interval [a, h] C Vt. By the definition of Vt and compactness, it is clear 
that there exists 5 = 6a,b, so that v' < —6 in V'i^a, h) sense. 

Fix a function z G C^{a,b). Clearly, for each such z, there exists Aq = \o{z,6a,b), so 
that for all < A < Aq, {v + Xz)/\\v + A2;||l9 satisfy all constraints. That is, we do not 
need to restrict over z in terms of positivity or bell-shapednes^. Thus, we have that for 
all < A < Ao(^), 



Je 



v + Xz \ _ Je{v + Xz) _ J^^'' + 2X{Mv,z) + 0{X 



\v + Xz\\l,J \\v + Xz\\l, {1 + Xqf^^vi-^{x)z{x)dx + 0{X^)y/i 

jmax ^ 2X{MV, Z) + 0(A2) 

1 + 2A /_!^ v''~^{x)z{x)dx + 0(A2) ~ 
= + 2X{{Mv, z) - z)) + O(A') 

Since ( \\J+xlu, ) ^ Jf"^^ conclude that for all z e C^{a, h) 

(19) (A^t; - J,°^""t;'-\ 2) < 0. 

Since there are no rectrictions on z (other than compact support), it follows that 
{Mv — Jf'^^v'^~^, z) = and hence 

(20) Mv - J^'^^^t;"-^ = 0, 

in the V{a,b) sense. This is the Euler-Lagrange equation that we were looking for. Note 
that according to our derivation, it holds in the compact subsets of Q. Note that from 

here, we obtain v = (^-j^Aiv^ ' ^ , which is smooth. These could be iterated further to 
show that V G C°° on open sets (a, b), whenever [a, b] C Q. 

3.2.2. Euler-Lagrange over non-trivial intervals ofOf". Our goal will be to show that such 
non-trivial intervals do not exist . To that end, we will first assume that they do exist and 
then, we will be able to derive an Euler-Lagrange equation on them, which will then lead 
to a contradiction. Take such an interval, say [aojC^i] C fi'^, < ao < cti. 

Claim: v = const, in V[aQ,ai] sense, i.e. v is constant on connected component^ of 

Proof. To prove that, assume the opposite, namely that v is not a constant in ©'(a^ai). 
Thus (since we know v' < 0), there exists a test function < ipo ^ "^(00,01), so that 
{v',ipo) = —a < 0. For some small k G (0,1), the set = {x G (ao,ai) : ipoi^) > f^} 
will be nonempty and open. Take an interval (fig, Si) C V^, so that ai — < 1- Take an 
arbitrary function ip G P(ao,ai), so that < ip{y) < k for ?/ G (ao,ai). It follows that 



"'which will be the case in the Section r3 . 2 .21 below 

^But it may be a different constant on a different component 
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ip{y) < K < ipo{y) and hence {v',ip) < {v',ipo) = —a < 0. Also J 'ip{y)dy < K{ai — Oq) < 1, 
whence 

(21) {y'^^)<-a<-a I ^{y)dy, 



for all ip G 'D'{aQ, ai), with {] < < n. One can now extend fl2Tl) to hold for all < ^/^ G 
P'(ao,ai). Hence, (ao,ai) C in contradiction with (ao,ai) C (ao,ai) C ^2^^. □ 

Now that we have established that v is a.e. constant on any non-trivial interval (ao, cti) C 
Vf^, we will derive the Euler-Lagrange equation for ( !T6|) on it. We will consider first the 
case ao > 0, the other case will be considered separately. 
Case I: ao > 

To fix the ideas, we consider first the case when the interval [ao, ai] is isolated from the 
rest of fi'^, that is, there exists r > 0, so that (ao — r, ao) C VL and (ai, ai + r) C VL. 

Fix 6^ > be so small that (ao — ^,ao),(ai,ai + 0) <Z Vl. Consider a test function 
z e Co~(ao-^,ai + 0). 

Clearly, since v = const on (ao, ai), we need to require that the function 2;' < on (ao, ai), 
in order for v + Xz to be bell-shaped functiorH (recall A > 0). Denote 

(22) G{x) =Mv- J^''''v''-\ 

Since v = const, a.e. on (ao,ai), G is a continuous function on (ao,ai). Following the 
approach of the previous section, we derive the equation f lT9|) . which in the new notation 
says that 

(23) {G,z)<0, 

for all z, so that v + Xz is an admissible entry for ( [TB]) . We will show that (12^ imphes 

(24) [ G{x)dx<0 



ai 



(25) / G{x)dx > 



for all b G (ao,ai). Applying for b = ai— and (!25|) for b = ao+ and taking limits, 
yields f^^ G{x)dx = 0. 



Claim: From and ( 125|) . one can infer that G = on (ao, ai). 

Proof. By the definition of G, the fact that v = const, and Lemma |3l we conclude that G 
is non- increasing and continuous function on (ao, ai). Now, we have several cases. If there 
exists a 6 G (ao, ai), so that J^^^ G{x)dx < and J^^^ G{x)dx > 0, then by continuity, there 
will be 60 £ {(iQ,b), so that G{bo) < and bi G (&, ai), so that G{bi) > 0, a contradiction 
with the fact that G is non-increasing. Otherwise, for all b G (ao,ai), we have that either 
/^^ G{x)dx = or /^"^ G{x)dx = 0. But since J^^ G{x)dx = 0, 

ai pai pb 

G{x)dx = / G{x)dx — / G{x)dx = — G{x)dx. 

J ao J ao J ao 



^so that it is acceptable entry for the minimization problem (jl6p 
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Thus, JI^^^ G{x)dx = for all b G (flo; cti)- Thus, for any Bq, hi G (flo; cli)-, we get ^^^^ G{x)dx = 
J^^^ G{x)dx — G{x)dx = and hence G{x) = in (oq, oi). □ 

Thus, the Euler-Lagrange equation is in the form G{x) = for x G {a^, Oi) or 

(26) Mv - J^'^'^v"-^ = 0, 

if we can show ( l24l) and ( l25l) . 

We show only ( l24l) . the proof of ( l25l) is similar. Fix b : aQ < b < ai. For all < 5 << 1, 
introduce an even test function, which for x > is given by 

( 0<a;<ao-25 
zs{x) = < 1 ao — 5<x<b + 5 
[ x>b + 26 

so that zs G C°^, zs is strictly increasing in (ao — 25, ao — S) and zs is strictly decreasing in 
(6 + 5,6 + 25). 

Note that v + Az^ is acceptable for the maximization problem ([T6|l foi@ < A < A5. 
According to f l23|) . {G,zs) < 0. Hence 

G{x)dx = limsup / G{x)zs{x)dx < 0, 

'ao <5-).0 Jao-25 

which is (El). 

The general case, in which [ao, ai] is not isolated from (i.e. there is no r > 0, so that 
(oq — r, ao) U (ai, Oi + r) C Q) is treated in the same way. Indeed, this was needed only in 
the very last step, in the construction of the function zs- But clearly, one can carry out a 
similar construction of zs, if one has a sequence of intervals Is^ = {cs^ —Sj, cs^) C Q, cg. < ao, 
so that limc^^ = oq. Similarly, for the proof of fl25l) . one needs a sequence of intervals 
Jsj = {dsj,ds^ + Sj) C Q, ds. > ai, so that lim^^. = ai. 

Finally, it remains to observe that every non-trivial interval of is contained in [ao, ai] C 
Qf^ with the property dist{aQ,VL) = = dist{ai,VL) and hence, we can carry the construc- 
tions of Zs and hence the validity of fl2^ and fl25l) follows. We then derive fl26|) on every 
such interval (ao,ai). 
Case II: oq = 

This case is similar to the previous one. We again assume that there exists r > 0, so that 
(— ai — r, — fli) U (ai, ai + r) C Vt, the general case being reduced to this one by arguments 
similar to those in the case ao > 0. 

Again, G is even, we have that v = const, on the interval (— ai,ai) and G is non- 
increasing in (0, ai). We will show that for every 6 G (0, ai), we have 

(27) / G{x)dx > 0, 



Jb 

(28) [ G{x)dx < 



Here even though zs is increasing in (ao — 26, ao — 5), this is acceptable since (ao — 25, ao — (5) C and 
hence, we have no restrictions over the test functions, as long as << 1 
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Let us first prove that assuming fl27|) and fl28l) . one must have G = 0. Indeed, if we apply 
([57]) for 6 ^ 0+ and ^ for b ^ ai-, we see that G{x)dx = 0. 

Again, if we assume that /o G(x)cix = for all b G (0, ai), we again conclude by the 
continuity of G that G{x) = : < x < ai. If one has Jg" G{x)dx < for some 
bo e (0, Ci), we have that /^"^ G{x)dx > (since J^^ G{x)dx = 0) and hence, G cannot be 
non- increasing function (as in Case I), a contradiction. Thus G{x) = in (0,ai). 

Thus, it remains to show ([27]) and (EH]). Fix 6 > and let < (5 << 1. For ([2H]), select 

( X <-b-26 
zs{x) = < 1 -b-6<x<b + S 
[ x>b + 26 

and zs is C°°, even and increasing in (—6 — 25, —b — 6) and decreasing in [b + 6,b + 26). 
Clearly, for A > 0, f + Xzs is admissible for f lT6]) . for some small A > 0. Thus, applying 
f[23]) for this 2:5 and passing to appropriate limits as 5 — t- 0+, we obtain, as above 



2 / G{x)dx = / ^(x)^^ = liminf / G{x)zs(x)dx < 0. 
Jo J-b '5-^0+ J 

For ( 127]) . we construct 2;^ to be an even function, so that 

( X <b-26 
zs{x) = < —1 b — 6<x<ai + 5 
[ x> ai + 26 

Now, we require that zs is decreasing in {b—26, b—6) and it is increasing from {ai+6, ai+26). 
Note that zg is still acceptable as a perturbation - in the sense that v + A^^ is non- increasing 
in (0, 00) for all A = A(5) > small enough (for this, recall that (ai + 6,ai + 26) C Q and 
hence, there exists <j{S), so that v' < —a in V'^ai + 6,ai + 26) sense). We have again 

— / G{x)dx = limsup / G{x)zs{x)dx < 0. 

Jb <5-i'0+ J 

Thus, we have established f lTT]) and f[25]) and thus the Euler-Lagrange equation 

(29) = G{x) =Mv- J^''\!''-\ 

3.2.3. The set O,'^ consists of isolated points only. We will now show that if an equation 
like ( |29l) holds in a non-trivial interval, say (a, b) and f is a bell-shaped, locally integrable 
function, which is constant on (a, 6), then v{x) = const, a.e. on (which would be a 
contradiction). 

Indeed, A4v is in fact a differentiable function on {a,b), which is non- increasing on 
(0, 00), according to Lemma |3l Thus, taking a derivative of f l29]) (and taking into account 
that V = const on (a, b)) leads to 

px+l px 

(30) = (Mv)'{x)= v{y)dy- v{y)dy. 

J X J x~l 

for all X G (a, b). 

If 6 > 1, we see that since v is non-increasing in (0, 00) (by Lemma [T]), it follows that 
X — )■ J^~^^ v{y)dy is non- increasing and continuous. By fl5CT]) . it follows that J^^^ v{y)dy = 
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const, for x G {a,b). Hence, we must have v = const, a.e. in every interval in the form 
(a; — 1, X + 1) (for all x G (a, b)). By iterating this argument in all (0, oo), a contradiction. 
If {a,b) C (0, 1), we can again argue as in Lemma |3] to establish that again v = const, in 
(0, oo). Thus, we cannot have non-trivial intervals {a,b) C Q'^ and hence Q'^ consists of 
isolated points only. 

3.2.4. The Euler-Lagrange equation on R^. Before deriving the Euler-Lagrange equation 
for the maximizer v, let us recapitulate what we have shown so far for v. We managed to 
show that is a dense open set, so that fl'^ consists of isolated points only. Finally, on Q, 
f is a continuous function, and the equation fl20|) holds on every interval (a, b) C Q. 

We will now show that (l2Up holds for almost all x G Q'^. First of all, recall that 
G = A4v — j™'^'^'?-! is locally integrable function (as a sum of an L"? function and L~ 
functions) and hence, almost all points are Lebesgue points for it. 

Let Xo > 0, Xo G fi'^, so that Xq is a Lebesgue point for G. We have shown that Xq is an 
isolated point of Q^, which implies the existence of intervals inside Q, which approximate Xq. 
That is, there are {aj — 6j, aj + 6j) C {bj — 6j, bj + Sj) C fl~^, so that aj + 6j < xq < bj — 6j 
and limj aj = limj bj = xq. In addition, we can clearly select these intervals to be very short, 
namely we require limj ^-3^ = 0. Construct now a sequence of even C°° test functions, 
given by 

( <x <aj- 6j/2 
Zj{x) = < 1 ttj < X < bj 
[ x>bj+ 6j/2 

where zj is strictly increasing in (oj —6j/2, aj) C Q and strictly decreasing in {bj, bj + Sj/2). 
We have already shown that functions of the form v ± Xzj will be non- increasing in (0, 00) 
and it will otherwise satisfy all the restrictions of the optimization problem f lT6|) . provided 
< A << 1. Thus, accoridng to (fT^ . we have {G, zj) < and —{G, Zj) < and hence 



G{x)dx+ / G{x)zj{x)dx+ / G{x)zj{x)dx = {G, zj) = 

Jaj—Sj/2 Jbj 

Dividing both sides by bj — aj, and taking lim as j — 00 (noting that Xq is a Lebesgue 
point), we get 

(31) G(xo) + lim — ^^[^ f ' Gix)zAx)dx + ^ /" Gix)zJx)dx] = 0. 



Thus, we need to show that the limit above exists and it is equal to zero. To that end, 
note that H^jHoo < 1 and 



Jaj-5,/2 



G{x)zj{x)dx + — / G{x)zj{x)dx\ < \\G\\l^ < {\\v\\lcc + Jl 



max /riNq— 1^ 



< (t;(0) + ""t;(0)^ 



maxll ||<J-1 

e II^IIl°° 



since v is non- increasing a.e. in (0, 00). The last inequality, combined with lim, , = 
shows that (13T|) implies G{xq) = 0. Hence, for all Lebesgue points of G, we have G{x) = 0. 
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Thus, 

(32) Mv - J^^'^v'^-^ = a.e. 

3.3. Taking a limit as £ — )■ 0: Constructing a solution to (j5]). From (132|) . 

(33) Mve - J^^'^'vl-^ = 0, 

which is satisfied a.e. and also in T>'{0, e^^ — l) sense. From it, we learn that Vg is H^{0, e"^ — 
1) and consequently, by iterating this argument, H°°{0,e~^ — 1) = fl^i -^'^(O^ ^"^ ~ 1) 
function. Recall also that by construction ||fe||L<7(Ri) = 1 and there exists Eq > 0, so that 
inf£<£Q J^^"^ > 0, see f|T7|) . We will now take several consecutive subsequences of e — )■ 0, in 
order to ensure that the limit satisfies ([5]). 

First, take — )■ 0+, so that lim„ J™^"" = limsup^^Q J™^"" := Jq > 0. Second, out of this 
constructed sequence £„, take a subsequence, say En,., so that Ve„^ ^ t> in a weak L'^ sense, 
for some v > 0,v G L'^. This is possible, by the sequential compactness of the unit ball 
in the weak L"^ topolog}!^. By the uniqueness of weak limits (by eventually taking further 
subsequence), we also get v^~^ V^'^ in weak L~ sense. We also have A^fe„ ^ A4v in 
the weak L'^ topology, since for every test function ijj G C^, we have by the self-adjointness 
of M and M : U' U', 1 < r < oo. 

Thirdly, we show that the limiting function v is non-zero. To that end, it will suffice to 
establish that 

(34) liminf \\Ve\\Li(-e-Ke-^) > 0. 

Assuming that flMl) is false, we will reach a contradiction. Indeed, let 6j — t- 0+ be a 
sequence so that lim,- lU,/ r-i g-i) = 0. Thus, 



Jx-1/2 

1 7 / 



3 



We now use a refined version of the Gagliardo-Nirenberg estimate that we have used before. 

1/ vs^{y)dy\^dx = \ vs.^{y)X(-5-^-i 5-^+i)dy\^dx < 

-<5-i Jx-l/2 •'-^J^ Jx-l/2 ^ ^ 

< \\Q[v5,X(-5-'-1,S-'+1)]\\l^ 

and hence 

I|QKX(-5-1-1,5-i + i)]||l2 < \\Q[v5,Xi-S;^^l,5-'+l)]\\XV'^^\\Qh,Xi^^^^ 

While a simple differentiation shows that on one hand, 

||<5KX(_5-i-i,5-i + l)]IUi.9 < ^vs^h^ = 2, 



But again, this argument, so far, does not guarantee that v ^ 
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we also have by Cauchy-Schwartz 

IIQKX(-5-i-i,5-i+i)]lli? < j {j V5,iy)X(-5~^-i,5;'+i)iy)dyydx < 



oo Jx-1/2 



5^ -1 J ^ 



The last inequality here follows by the fact that vs^ is non-increasing in (0, oo) and therefore 
have proved that 



J li Vs {y)dy < /o^ 'Vs.{y)dy, if > 1, which we have assumed anyway. Thus, we will 



liminf Jr^" < 0, 



which is in contradiction with (fTTll . Thus, we have established ( !34l) . 

We are now ready to take a limit as e — )■ in (155]) . Indeed, take (1551) for e = Sn^, k = 
1,2,.. .. Fix a test function ip e C^. There exists ko, so that for k > k^, supp ip C (0, Sn^)- 
Thus, for k > k^, we gelH 

Take a limit as A; — t- oo. By our constructions, we have that J^^^{v^~^ — Jo{v'^~'^,iIj), 

since J^^J" — )■ Jo, {v^^^,i/j) — )■ {v^^^^ip) by the weak convergence. 

We also have {J^Vs„^,ip) = {ve„^,J^4') -^k {v.J^ip) = {Aiv,ip). Thus, we have estab- 
lished the desired identity 

(35) Mv - Jov''-^ = 

valid for all x > 0. By the symmetry, it is also valid for x < 0. It is clear that v is now 
infinitely smooth^ function on (0, oo). Recall though, that for Theorem [1] we needed to 

solve JHw — w'^^^ = 0. One can easily construct w, based on the solution v of (!35|) . More 

]_ 

precisely, if we take w := Jq ^~''v, then w will satisfy ( |T5|) . Theorem [T]is proved. 

4. An ALTERNATIVE PROOF OF THE FRIESECKE-WATTIS THEOREM 

We quickly indicate how our ideas can be turned into a new proof of the Friesecke- 
Wattis theorem. Specifically, as we saw in the previous section, it is clear that if we are 
just interested in the existence of traveling waves for (but not in bell-shaped solutions), 
it is a good idea to consider following constrained maximization problem (compare to (fT6|) ) 



(36) 



Jeiv) = /f^-i \Qv{x)\'^dx — )■ max 
subject to V > 0; I{y) = f^i v'^{x)dx = 1. 



First off, the arguments in Section lSTTl apply unchanged (by just skipping the bell-shapedness 
of v) to prove that (l36l) has a maximizer, say v. 



%y testing ([33)) , which is vahd on the support of 4' 

'^starting with v G L''(R^), it is easy to conclude that ^Av is smooth, which in turn imphes that 
V G C\R^) etc. 
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Following the argument of Section 13.21 and more specifically, the following identities 

\\v + Xz\\% = \\v\\% + \q{v^-\ z) + O(A') 
J,{v + \z) = J,{v) + 2\{Q^v, z) + 0(A2), 

Je f II "IVh ) = -^^^^ + '^^{Mv, z) - J--(t;«-\ z)) + 0{X') < Mv), 

\\\v + Xz\\li J 

which were established there, it follows that 

{Mv - J'^^^'v'^-^ , z) < 
for all test function^ z. It follows that v satisfies 

Mv - J^^^'v'^-^ = 0. 

This of course produces a family Vs '■ supp C — + 1), which easily can be 

shown to converge0 (after an eventual subsequence) to a f : supp f C R\ which solves 

Mv - Jov''-^ = 0. 

1 

Setting w := Jq V again provides a solution to w^^^ = A4w as is required by ( IT5l) . 
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